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THE COMPRESSIBLE TURBULENT BOUNDARY LAYER ON A BLUNT 
SWEPT SLAB WITH AND WITHOUT LEADING-EDGE BLOWING* 

By James L. Hunt, Dennis M. Bushnell, 
and Ivan E. Beckwith 
Langley Research Center 

SUMMARY 

A finite-difference method has been developed to solve the equations for compres- 
sible turbulent boundary layers on swept infinite cylinders. Predictions by this method 
have been compared with experimental data on a 60° swept, blunt, slab configuration with 
and without leading-edge blowing. The test conditions were a free-stream Mach number 
of 8, a wall-to-total temperature ratio of 0.4, and a range of free-stream Reynolds num- 
ber based on leading-edge diameter Njj e>00 from 0.92 x 10^ to 9.3 x 10®. 

Three different approaches to the formulation of eddy-viscosity models for three- 
dimensional boundary layers were considered. Two of these formulations were used in 
the numerical solutions. Nonsimilar solutions for laminar flow were obtained by setting 
the eddy viscosity equal to zero. 

Predictions for heat transfer, surface streamlines, Mach number profiles, and 
boundary-layer thickness have been compared with experimental data for zero blowing. 
These comparisons indicated that for Nj^ e ,oo = 2.6 x 10®, the leading-edge boundary 
layer was turbulent, but as the flow expanded around the leading edge, the measured prop- 
erties of the boundary layer on the slab portion of the model approached those predicted 
for a laminar boundary layer. This apparent laminarization of the initially turbulent 
boundary layer was probably caused by the large streamwise pressure gradients and the 
rapid expansion of the flow in the leading-edge region. At the higher test Reynolds num- 
ber of 9.3 X 10®, transition back to turbulent flow was observed on the downstream part 
of the slab. Both eddy- viscosity models gave satisfactory predictions of the turbulent 
flow, provided the mixing length in the outer part of the boundary layer was taken as 0.05 
of the boundary-layer thickness. This fraction is generally larger on flat-plate turbulent 
boundary layers where values up to 0.10 have been used. The smaller values of the ratio 
required to obtain agreement between theory and the present data suggest that the level of 
turbulence in the present three-dimensional boundary layers was small compared with 
flat -plate flows at high Reynolds numbers. 

*This paper is an extension of paper no. 19 in NASA SP-228, 1970. 



The hypothesis of low turbulence levels in the leading-edge boundary layers of these 
tests was also consistent with the apparent laminarization phenomenon since the value of 
a parameter that characterizes the onset of laminarization was an order of magnitude 
below the critical levels observed for laminarization in two-dimensional flow. The small 
value of this parameter suggests that the turbulent boundary layer on the leading edge 
was highly susceptible to laminarization as would be expected if the initial turbulence 
levels were small. 

The measured heating rates with blowing indicated that the laminarization-like 
behavior and the subsequent transition back to turbulent flow were essentially independent 
of leading-edge blowing for the smallest blowing rates used in this investigation. How- 
ever, when the blowing rate was increased by about 40 percent, the measured heating 
rates remained at the higher turbulent levels for about 2 leading-edge diameters farther 
downstream. 


INTRODUCTION 

The boundary layers on large lifting vehicles are predominantly turbulent and three 
dimensional in nature. For such vehicles, the accurate prediction of frictional drag is 
essential. As the flight Mach numbers are increased above approximately 3, reliable 
predictions of heat transfer also become important, especially on windward and leading- 
edge regions where active or passive cooling systems are required. 

The boundary-layer flows on the windward and leading-edge regions of these large 
lifting vehicles can often be classified as belonging to a particular type of three- 
dimensional flow designated as quasi-two-dimensional since the governing equations can 
be expressed in terms of only two independent variables. If transition from laminar flow 
has occurred far upstream of the region of interest, the turbulent boundary layers on 
sharp cones at small angles of attack and on infinite cylinders at yaw belong to this class 
of flows. These flows may exhibit significant three-dimensional effects since the velocity 
vector within the boundary layer may have large components normal to the external 
inviscid streamline direction. The boundary-layer velocity component in this direction 
is generally referred to as the cross flow. 

Several previous methods developed for the calculation of compressible boundary 
layers of this quasi-two-dimensional class, as well as more general three-dimensional 
flows, were integral methods (refs. 1 to 5). The assumption of small (or zero) cross 
flow within the boundary layer was used in references 3 to 5. As shown by the results 
of reference 5, the assumption of small cross flow is generally more valid for turbu- 
lent flows than for laminar ones. This assumption requires that the magnitude of the 
velocity normal to the inviscid streamlines and the derivatives in that direction be 
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small; thus, it follows that the component of the surface shear normal to the inviscid 
streamline is small and has a negligible effect on the streamwise flow which develops 
independently of the cross flow. (See ref. 6, for example.) The Reynolds stress in 
the streamwise direction is then independent of the cross flow, and two-dimensional 
models for the Reynolds stress are directly applicable. 

For large streamline curvature and large surface blowing, the cross flow in turbu- 
lent boundary layers may be large. With large cross flows, the shear stress normal to 
the inviscid streamlines influences the streamwise flow, and the two-dimensional models 
of the Reynolds stress may no longer be adequate. Thus, boundary- layer solutions in the 
presence of large cross flows necessitate the formulation of models for both components 
of Reynolds stress (streamwise and cross flow) that will properly account for the three- 
dimensional nature of the flow. A finite -difference procedure for solving the equations 
for incompressible turbulent boundary layers without the small-cross-flow assumption 
has been developed by Nash (ref. 7) who assumed that the magnitude of the shear stress 
is proportional to the turbulent intensity. 

A finite-difference method is utilized herein which is applicable to compressible 
boundary- layer flows in the presence of large cross flows and surface blowing. This 
method of solution consists of the extension of the finite -difference procedure of refer- 
ence 8 to include the spanwise momentum equation for swept infinite cylinders. The com- 
putational advantages of two independent variables are thereby retained, while the bound- 
ary layer may have large cross-flow components. The computations are not carried out 
in streamline coordinates, and the only restrictions are those arising from the. assumed 
models for the components of the Reynolds stress in the chordwise and spanwise direc- 
tions and for the turbulent flux of enthalpy. 

The purpose of this investigation is to evaluate the proposed Reynolds stress models 
by comparison of theoretical predictions with new experimental data obtained on a 60° 
swept slab with and without homogeneous blowing at the leading edge. The test condi- 
tions were a free-stream Mach number of 8, a wall-to-total temperature ratio of 0.4, and 
free-stream Reynolds numbers, based on leading-edge diameter, from 0.92 x 10® to 
9.3 x 10®. Several solutions for test Reynolds numbers of 2.6 x 10® and 9.3 x 10® have 
been obtained with various combinations of two eddy-viscosity models and different 
mixing-length functions. The results are compared with measurements of wall heat 
transfer, oil-flow lines (representing surface streamlines), Mach number profiles, and 
boundary-layer thickness. The region of special interest is just downstream of the blunt 
leading edge where the boundary-layer profiles are far from equilibrium because of the 
upstream conditions of blowing and large cross-flow pressure gradients. 

This paper extends the results and analysis of reference 9 to include additional data 
at higher blowing rates for a Reynolds number of 2.6 X 10®. 
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SYMBOLS 


area 


A function of blowing parameter (see ref. 8) 

a i,bi, c i,di, e i,fi temperature coefficients in heat -balance matrix (34) 

(pv). 


B 


normalized blowing rate, 


'w 


(pu)< 


Cf 


• Ttjf 

skin-friction coefficient, •= =■ 

fPe u e 

specific heat of model material 


specific heat at constant pressure 


d 

G 


leading-edge diameter 


eddy-viscosity function defined by equation (19) 


H 

h 


constants in heat-balance matrix (34) 


total enthalpy, h + 


static enthalpy 


u^ + w^ 


heat-transfer coefficient 


k 

L 

l 

M 


molecular thermal conductivity 


reference length 


mixing length (eqs. (8) and (9)) 


Mach number (based on local velocity vector q) 


m 


mass flow rate per unit area 


4 


I 


N 


Pr,T 


N 


Re,s 


N 


Re,°° 


N 


St 


n 

P 

q 

q 

R 


T 

t 

u,v,w 


c n e 

turbulent Prandtl number based on static enthalpy, -tL- 

K 


reference Reynolds number, 


P sV 2H e L 
Ms 


free-stream Reynolds number based on leading-edge diameter, 

Stanton number, 77- — ■ — rv 

Pe u e( h aw “ h ) 

exponent in transformation (eqs. (26)) 
pressure 

heat -transfer rate per unit area per unit time 
local velocity vector, + w 2 
ideal gas constant 

distance along an inviscid streamline 

absolute temperature 

time 

velocity components in x-, y-, and z-directions (fig. 1) 


PooMoo d 

Par, 


V = v + ol 

p 


x,y,z 

a 

AY 

Ay 


volume 

Cartesian coordinate system (fig. 1), also coordinate system used in heat- 
conduction problem (fig. 21) 

2m AY(pRT) w 


2 . 


permeability, ft 2 , 

Pc r w 

wall thickness of porous leading edge 

wall thickness at thermocouple stations on slab 
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M 

6 boundary -layer thickness, usually evaluated at the point where ~ 0.995 

M e 

e eddy viscosity 

e* emissivity 

6 angle between z -coordinate direction and local velocity vector 

k eddy thermal conductivity (eq. (7)) 

A leading-edge sweep 

jit molecular viscosity 

t] transformed coordinate system (eqs. (26)) 

p mass density 

a Boltzmann constant 

r shear stress 

Subscripts: 

aw adiabatic wall 

c coolant conditions inside porous leading edge before injection 

e local "edge" of boundary layer 

g gas 

m model material 

max maximum 

min minimum 

n direction parallel to surface and normal to local inviscid streamline 
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o initial condition at time zero 

p direction parallel to the surface and parallel to the local inviscid streamline 

r local reference condition, may be either wall or edge 

s stagnation line 

T turbulent 

t total 

w wall or outside surface of model 

x,y,z in the x-, y-, or z-direction 

°° free-stream flow in tunnel 

A bar over a symbol indicates a mean quantity. A prime indicates a fluctuating 
quantity. 


ANALYTICAL METHOD 

Numerical solutions to the partial differential equations for the boundary layer on 
swept infinite cylinders are obtained by an implicit finite -difference procedure. The 
method of reference 8 for two-dimensional and axisymmetric flows is extended to swept 
cylinders by including the spanwise momentum equation. The Reynolds stress terms in 
the momentum equations are formulated by the use of eddy-viscosity concepts. The 
turbulent flux of static enthalpy is related to the Reynolds stress by the turbulent Prandtl 
number. 

The formulation of models for the eddy viscosity in three-dimensional boundary 
layers is of primary concern herein. Two different models for the eddy viscosity are 
utilized in the calculations: (1) the "invariant turbulence" model, in which the total shear 
is treated as a vector and the eddy viscosity is assumed independent of direction and 
(2) Lettau's vorticity-transfer hypothesis (ref. 10), in which the change of vorticity fol- 
lowing a turbulent eddy is expressed in a general vector form. 
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Boundary- Layer Equations 

The coordinate system for the idealized flow model is shown in figure 1(a). The 
instantaneous flow variables in the general conservation equations are divided into mean 
and fluctuating parts, and the Reynolds averaging process is applied. The following equa- 
tions are then obtained after the higher order terms are neglected according to the usual 
order-of-magnitude analysis for thin boundary layers: 

Continuity 

^pu) + ±(py) = 0 (1) 

Chordwise momentum 


— 9 a — 9u 9p d (- da -~TT\ 

pu + pv -ft = + — p. -ft - pv'u' 

ax H ay ax ay ^ ay / 

Spanwise momentum 

— cfw 9w a f— 9w — : — A 

pu ^ + pv ^ = ^i M ^r- pvw 7 


( 2 ) 

(3) 


Energy 


pu® + pv*U± 

ax ay ay 


^ - p v'h' - puv’u’ 

C P ay 


— — t — 7 - a 

pwv’w’ + u — 

9y 


- 2 — 2 N 

r u + w 


(4) 


These equations are valid in the entire stagnation region of a cylinder if a characteristic 
Reynolds number, based on the radius of curvature of the leading edge, is large. (See 
ref. 11.) 

Before equations (1) to (4) can be solved, expressions must be supplied for the 
Reynolds stress or turbulent shear terms in the momentum equations and the turbulent 
flux of enthalpy in the total-enthalpy equations. The approach used in the present method 
is to model these terms as functions of the mean-flow variables. 


Eddy Diffusivity Coefficients 

The concept that the Reynolds stress in turbulent flow is proportional to a momen- 
tum exchange coefficient times the mean-flow velocity gradient was first proposed by 
Boussinesq, as discussed in reference 12. This concept is based on an assumed analogy 
between eddy viscosity and molecular viscosity. The shear components in the chordwise 
and spanwise directions are then written as 
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SWEPT INFINITE CYLINDER 


w e =w »:^ 


( )=o 




(Id) Test model configurations. 

Figure 1.- Basic model configuration and coordinate definitions. 
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9y 
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— 9u 
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9y 
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where the eddy viscosities in the x- and z-directions might in general be different. Since 
the total shear can be treated approximately as a vector quantity in boundary -layer flows, 
its magnitude may be written as 



The eddy thermal conductivity k is defined in a manner similar to that used in 
equations (5) for the eddy viscosity. The total heat flux in equation (4) is then expressed 
as 


q = — - p v'h’ 

c p ay 


k ah [ k ah 

c p 9y Cp 9y 


(7) 


The remaining two fluctuating terms in equation (4) are products of velocity and Reynolds 
stress; thus, they are modeled according to equations (5). 


Eddy- Viscosity Models 


The simplest approach to the formulation of models for the Reynolds stress is 
based on Prandtl’s mixing-length hypothesis (ref. 13). For two-dimensional flow, this 
hypothesis states that 


u- cc-z^i 

9y 


v’ «z 


9u 

9y 


( 8 ) 


The turbulent shear and eddy viscosity are then 


t t s -pv’u’ = pi 2 


du 

9y 


3u 

9y 


-2 

e = pi 


au 

9y 


( 9 ) 

( 10 ) 


where the quantity l is the mixing length, which is presumably related to some charac- 
teristic length of the turbulent motion. In the near-wall region of a boundary layer, l is 
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assumed to be proportional to the distance from the wall. In the far-wall region or in 
free turbulence, l is assumed to be proportional to the width of the mixing layer, that 
is, to the boundary-layer thickness or the width of the mixing jet. Although the details 
of such a transfer mechanism are not well understood, the basic concept gives satisfac- 
tory results even at hypersonic Mach numbers in the presence of large heat transfer and 
large pressure gradients. (See refs. 8 and 14.) One purpose of the present paper is to 
determine whether the basic mixing-length concept can be extended to three-dimensional 
flows. Three different approaches to this problem are considered. 

Independence principle .- The chordwise development of the laminar, incompressible 
boundary layer on a swept infinite cylinder is independent of the spanwise flow. (See 
ref. 15.) This concept is referred to as the "independence principle." The application 
of this principle to turbulent flow was first attempted by Young and Booth. (See ref. 16.) 
Although the variable density in compressible flow couples the chordwise and spanwise 
momentum equations, it is of interest to consider whether a valid mixing-length expres- 
sion for the Reynolds stress can be formulated from the independence principle. Thus, 
if the chordwise flow is independent of the spanwise flow, the chordwise shear would have 
to be of the form 


t x,T 



&u 9u 

9y 9y 


(ll) 


The simplest expression for the spanwise shear that retains the correct form at the stag- 
nation line is 


x,T 



8w 9w 
9y 3y 


( 12 ) 


The resulting expressions for the eddy viscosities are 




9u 

3y 

9w 

Qy 


(13) 


Hence the eddy viscosities are different in the two directions and therefore depend on 
coordinate direction. The turbulent Prandtl number would then presumably be based on 
the "total" eddy viscosity, with the result that 


N Pr,T 



(14) 
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However, substitution of equations (13) into equation (6) gives 



( 15 ) 


and it is evident that a total eddy viscosity like that of equation (14) cannot be obtained in 
explicit form from equation (15). The requirement of the independence principle that the 
eddy viscosity should depend on coordinate direction may therefore be erroneous. Exper- 
imental data of Ashkenas and Riddell (ref. 17) also indicate that the independence principle' 
may not apply to turbulent flows. (See also discussion in ref. 18.) In view of these dif- 
ficulties, the formulation of eddy viscosities from the independence principle will not be 
pursued further herein. 


Invariant turbulence .- An alternate approach to the independence principle is based 
on the concept that the eddy viscosity should depend only on the properties of the turbu- 
lence and a local eddy scale as in the methods of Prandtl (ref. 19) and Glushko (ref. 20). 
The application of this concept to a mixing-length approach suggests that the eddy vis- 
cosity would be a scalar function independent of the coordinate direction. Accordingly, 
the components of the turbulent shear terms are 


X >T ay ay 


_ -,2 9G aw 

t z,T = P l — — 




ay 9y 


(16) 


and the eddy viscosity is then 


-,2 8G 

e = e x = e z = pl — 


(17) 


z ay 

To determine the scalar function G, equation (17) is substituted into equation (6) to obtain 

_- x 2 ,_\jni/2 




: )T 


(18) 


Then by analogy with equations (9) and (10) , where the velocity-gradient function for the 
shear component is repeated in the eddy viscosity, the expression for G is 
~-\2 ,^ 1/2 

(19) 


8G. 

9y 


/ au\ /aw \ 

\dy) \dy)_ 


The final expression for the eddy viscosity then becomes 


2 / .. — \2 


e = pr 


(I) + (f ) 


|l/2 


( 20 ) 
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Bradshaw (ref. 21) has derived a set of differential equations for the two components 
of shear stress based on the Navier-Stokes equations and several assumptions regarding 
the behavior and formulation of correlation terms for fluctuating pressure and velocity. 
These equations indicate that the directions of the shear stress components are not gen- 
erally the same as the directions of the corresponding mean velocity components. How- 
ever, in the near-wall region, Bradshaw's equations reduce to the same form as the pres- 
ent results (eqs. (18) and (19)). 

Lettau's vor ticity -transfer hypothesis .- In reference 10, Lettau proposed a modifica- 
tion of Taylor's vorticity -transfer theory, which accounts for the combined action of con- 
servation and adjustment of the vorticity of a fluid element as it is convected along an 
eddy trajectory. The resulting expressions for the three components of fluctuating veloc- 
ity reduce to simple forms for the present application. That is, by the use of conventional 
order-of-magnitude analysis for a thin boundary layer on a swept infinite cylinder, these 
expressions for the fluctuating velocity components reduce to 



For simplicity, the various components of an "eddy displacement vector," represented in 
equations (21) by l x , ly, and l z , are assumed to be independent of direction. The com- 
ponents of the turbulent shear then become 



The eddy viscosity is then given by 

which is the same for both shear components, but in contrast with the invariant -turbulence 
model of equation (20), this model is the scalar sum of two vector quantities. Recently, 
Lettau (ref. 22) has generalized his original hypothesis of reference 10 in an attempt to 
account for possible variations of his eddy displacement vector and has applied the results 
to a free -turbulence problem. 
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Mixing-Length Function 


The mixing-length function used for the present calculations is 

l 

6 

This expression is the same as the function used in reference 8 except that the gas proper- 
ties in the exponential are herein always evaluated at the surface, and the function f(y/5) 
is modified to account for smaller levels of turbulence. That is, the fj function of ref- 
erence 8 (table I in ref. 8) was modified by reducing the maximum value of the function 
for y/6 > 0.1 from 0.09, to 0.07 or 0.05. The Prandtl slope in the near-wall region 

df ic 

— = 0.4 was retained. Also, the functional dependence of A on the blowing param- 

d(y/6) 

eter 2B/Cf as given in figure 1 of reference 8 has been used in the present calculations. 

Variation in Turbulent Prandtl Number 

The variation in turbulent Prandtl number through axisymmetric compressible 
boundary layers with large heat transfer and large favorable pressure gradients is c 
treated extensively in reference 8. The results indicated that values of N pr p = -E- 

(the "static" turbulent Prandtl number) in the outer part of these boundary layers may be 
as small as 0.5. Values of Np r «p obtained from experimental data on flat plates and 
reported in references 23 and 24 are generally greater than 1.0 near the wall and tend to 
decrease below 1.0, or in some cases to nearly 0.5, farther from the wall. On the basis 
of these results, the Np r <p variation used in all solutions reported herein was taken as 
a simple ramp function of y/6 as given by the following relations: 

(25) 

Some check calculations with the alternate assumption of Np r p = 0.9 reduced the pre- 
dicted heat transfer by at most 15 percent, which was considered to be within the range 
of experimental uncertainties. For the present conditions, the choice of the Np r p 
function is apparently not critical. 

Numerical Procedure 

The Reynolds stress terms in equations (2) and (3) are replaced by the eddy- 
viscosity expressions (eq. (20) or (23)) with the mixing-length function given by equa- 
tion (24). These eddy-viscosity and mixing-length expressions are also used in the 



p w y 

|1 - exp(- — 
A M-yy 



f(| 

o 


( 24 ) 
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energy equation (4) after replacement of the turbulent heat flux by the eddy-conductivity 
relation of equation (7) and the introduction of Np rj >p. The resulting equations are trans- 
formed to the coordinate system £,17 defined by 


5 


V 


(l) - N Re,s I 


x/L (pmL u 


0 (PM)s \j2Hl 


H) 





Re,s 


Ue/^H^ ry/L d 

( 2 |)n Jq p s 



(26) 


where n is adjusted to obtain a nearly constant boundary -layer thickness in the 
£,77 plane in order to increase the computational efficiency. 

The numerical procedure used to solve the transformed equations is an implicit 
finite-difference procedure similar to that of reference 25. The partial derivatives are 
replaced by linear -difference quotients. The result is a set of N - 1 linear algebraic 
equations for each of the unknown velocity components and enthalpy (u, w, and h) at the 
N grid points for the next downstream station. The matrix for each set of these linear 
equations is tridiagonal, so that an efficient algorithm is used in their solution. After 
the solution for u and w is obtained from the momentum equations, the transformed 
normal velocity is obtained from the continuity equation; then the enthalpy (at each grid 
point) is obtained by the same numerical method used for u and w. This procedure 
is repeated with updated values used in all difference coefficients until convergence is 
obtained. 

The solutions were started at a small but finite chordwise distance from the stag- 

_ . H - H w 

nation line with an assumed set of profiles for u/u e , w/w e , and — — ■ — . This pro- 

H e - H w 

cedure was used for the present solutions because the correct limiting forms of the equa- 
tions at x = 0 had not yet been incorporated in the computer program. It was found that 
consistent solutions for all variables could be obtained within a distance equivalent to 2 
or 3 boundary-layer thicknesses by using a large number of small steps in the chordwise 
variable £. These solutions depended on the particular eddy-viscosity model used but 
were independent of the input profile shapes and input boundary-layer thicknesses, pro- 
vided reasonable assumptions were used for these quantities. 

The implicit finite-difference procedure as just described can provide, for all 
practical purposes, nearly "exact" solutions to the partial differential equations since 
the method is numerically stable. Of course, the accuracy of the solutions depends on 
the step sizes, the convergence criteria for the iterative cycles, and the accurate speci- 
fication of boundary conditions. Essentially by trial-and-error processes, these matters 
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have been treated during the course of this and previous investigations (refs. 25 and 8) 
with the result that the present solutions are believed to be accurate to within 1 percent 
or better. Consequently, the comparison of theoretical predictions with data provides an 
indirect evaluation of the eddy- viscosity models, the mixing-length functions, and turbulent 
Prandtl number distributions since for a given set of boundary conditions, these are the 
only ingredients in the theory that can account for different predictions. Thus, by using 
various eddy-viscosity models and mixing-length functions in the theory and comparing 
results over a range of Reynolds number and with as many different experimental mea- 
surements as possible, the optimum formulation of Reynolds stress and turbulent Prandtl 
number can be determined. Such a process may be considered as "numerical experi- 
mentation” to determine the best form and values of the Reynolds stress. 

EXPERIMENTAL METHODS 

Experimental data have been obtained on a 60° swept slab with and without blowing 
at the blunt leading edge. The surface heat transfer, Mach number profiles, boundary- 
layer thickness, and surface streamlines were obtained for the model without blowing. 

For the model with blowing at the leading edge, the only data available are for surface 
heat transfer on the leading edge and downstream on the slab. The wind-tunnel facility 
and models will be described in detail in the following sections. 

Facility 

All the data presented herein were obtained in the Langley Mach 8 variable -density 
hypersonic tunnel. This blowdown-type facility has a contoured axisymmetric nozzle with 
an 18 -inch-diameter test section. Transient testing techniques, that is, rapid exposure 
of models to established flow conditions, are accomplished by means of a model injection 
mechanism located directly beneath the test section. Windows are located on both sides 
and at the top of the test section for lighting and photographing models. The free-stream 
Mach number varies slightly with free-stream Reynolds number (ref. 26). Data will be 
presented herein for free-stream Reynolds numbers, based on leading-edge diameter, of 
0.92 x 10^, 2.6 x 10^, and 9.3 X 105. The corresponding free-stream Mach numbers used 
in the data reduction are 7.81, 7.94, and 7.98. The stagnation temperatures are 1285° R, 
1390° R, and 1475° R, respectively. 

Models and Instrumentation 

For this investigation the test configuration is a 60° swept slab with a leading-edge 
radius of 1/2 inch. One pressure model, three plastic heat -transfer models, and one 
porous-leading-edge model were used in the experimental tests. A sketch of the basic 
configuration used for the pressure model is shown in figure 1(b). The coordinate system 
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used to locate the instrumentation on the various models is also given in this figure. The 
forward portion of all models (ahead of the dashed lines) was geometrically identical with 
the pressure model shown. The aft portion of the plastic models and the porous-leading- 
edge model were modified as indicated by the dashed lines in the figure. The mounting 
plate as shown in the sketch was used for the pressure model and the plastic heat -transfer 
models. The mounting plate for the porous -leading-edge model was modified, as will be 
shown subsequently; however, the leading edge of the plate was 1/4 inch ahead of the apex 
of all models. 

Pressure model .- The pressure model was made of stainless steel and was pro- 
vided with seven chordwise rows of pressure orifices over a 12-inch span as shown in 
figure 2. The distribution of these orifices in the three chordwise rows used in this 
investigation is given in table I. A boundary-layer pitot -tube rake is shown attached to 
the model in figure 2. Details of the rakes used will be given in a subsequent section. 

Phase-c hange heat-transfer models .- A photograph of one of the three plastic 
models used with the phase-change coating technique is shown in figure 3. The model 
shown was cast from an epoxy resin with a silica base. One of the other plastic heat- 
transfer models was cast from a slightly different mixture of epoxy and silica, and for 
the third model, mica was used as the base material. Therefore, these three models had 
different thermophysical properties. The value of the combination of these properties 
required in the phase-change method is \jpc k. The values of this quantity for the three 
models were 0.073, 0.075, and 0.034 Btu/ft^ °R sec*/^, respectively. 
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Figure 3.- Plastic model used with phase-change technique. 


L-71-512 


Porous-leading-edge m odel.- A photograph of the porous-leading-edge model used 

in this investigation is shown in figure 4. The leading edge of this model consisted of a 

porous hemicylinder 10- inches long and approximately 3/32 inch thick. The hemicylin- 

0 

der was made of sintered stainless steel with a porosity of approximately 13 percent. 

The leading edge was electron-beam welded to a swept slab which consisted of a 
machined, ribbed structure with a thin skin. The mounting-plate configuration for this 
model is different from that shown in figures 1(b), 2, and 3. In particular, the forward 
section of the plate was machined into an elliptical planform which projected beyond the 
periphery of the model apex by 1/4 inch. 

The coolant air was fed into a chamber just behind the porous leading edge. The 
chamber could be pressurized to 415 psia. The pressure and temperature of the coolant 
air were measured in this chamber. 

Thermocouples were installed in the porous leading edge and attached to the thin 
skin of the slab. The locations of the thermocouples used in this investigation are given 
in table II. To install the thermocouples on the porous leading edge, 0.0 16 -inch-diameter 
balls were formed at the thermocouple junctions. Teflon-coated chromel-alumel wires 
of 0.003-inch diameter were used. The thermocouple wires were then inserted into a 
0.013-inch-diameter hole which was burned by an electric-arc technique through the 
porous leading edge at an angle of 45° to the surface. The axes of the holes were alined 
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L- 71- 513 

Figure 4.- Porous-leading- edge model. 


in spanwise planes which in turn were normal to the surface. The thermocouple junction 
(a 0.016-inch-diameter ball) was tapped into the 0.013-inch-diameter hole and finished 
flush with the outer surface. The hole was then filled with epoxy resin to prevent direct 
leaks of coolant air around the junction. This procedure resulted in a thermocouple junc- 
tion within 0.004 inch of the external surface. Good thermal contact of the junction with 
the porous matrix was obtained. 

On the slab portion of the model the thermocouples were welded to the back surface 
of the stainless-steel skin, which was 0.030 inch thick. 

Boundary-layer pitot -tube rakes .- One of the boundary- layer pitot -tube rakes used 
to obtain Mach number profiles is shown attached to the pressure model in figure 2. A 
photograph of these rakes showing the tube assembly is presented in figure 5(a) . The 
rakes consisted of two rows of three pitot tubes each. Sketches of the pitot -tube con- 
figurations with pertinent dimensions are given in figure 5(b). 

Test Techniques and Data Reduction 

All data were obtained by using a transient testing technique in which the tunnel 
was started and brought to the desired operating conditions and then the model was rapidly 
injected into the airstream by a pneumatic piston. An analog-to -digital data recording 
system with a maximum resolution of 40 points per second was used for thermocouple and 
static -pressure data. 
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L-69- 7083.1 

(a) Photograph of raises . 


SIZE OF ORIGINAL CIRCULAR TUBING 

I.D.,in. O.D., in. 
SMALL RAKE 0.030 0.050 

LARGE RAKE .040 .060 



SMALL RAKE LARGE RAKE 


(b) Sketches of pitot tubes. Dimensions are in inches. 
Figure 5-- Pitot pressure rakes used, for boundary- layer surveys. 
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Pressure data .- Strain-gage pressure transducers were used to measure the static 
pressures on the model shown in figure 2. The locations of the pressure orifices (0.040- 
inch diameter on the leading edge and 0.060-inch diameter on the slab) used in this inves- 
tigation are given in table I. The output of the transducers was monitored to determine 
when time -steady pressures were reached, then the data were recorded on the readout 
system at the rate of 20 points per second. 

The pitot -tube rakes shown in figure 5 were mounted on the pressure model as 
indicated in figure 2. Pitot-pressure data were obtained at five chordwise locations on 
the slab portion of the model. The rakes were alined so that the pitot tubes were parallel 
to the free-stream flow direction. The plane of the pitot orifices was located approxi- 
mately 1/4 inch behind a static-pressure orifice at each chordwise location. The 
y-location of the tube assembly was adjusted with shims between the rake and model. 

The pitot -pressure data were read from photographs of a mercury manometer board. 

The ratios of the local surface static pressure to the pitot pressure were used to 
compute the Mach number distributions through the boundary layer. The boundary- layer 
thickness was taken as the distance from the surface where the Mach number became 
approximately constant. 

Oil flow .- In order to obtain the direction of the wall streamlines and the overall 
flow pattern, small dots of a mixture of oil and lampblack were placed on the slab portion 
of the pressure model. The model was then rapidly injected into the airstream. After 
exposure to the flow of 3 to 8 seconds, depending on the Reynolds number, the model was 
retracted and the resultant patterns photographed. 

Heat transfer for no blowing .- Heat -transfer-coefficient distributions on the plastic 
"stycast" models (fig. 3) were obtained with the phase-change coating technique (ref. 27). 

In this technique, a thin surface coating of material that undergoes a visible phase change 
at a known temperature is sprayed on the model. The model is then injected rapidly into 
the established flow and the phase-change patterns which develop as the model is heated 
are recorded by motion-picture photography. The patterns of isothermal lines so obtained 
may be transformed to lines of constant heat -transfer coefficient provided the distribution 
of adiabatic wall temperature is available and the thermophysical properties (pck) of the 
model material are known. Exact solutions of the heat-conduction equation for the specific 
geometry are not generally practical; therefore, the local heat -transfer coefficients for 
the model are determined from the solution for a semi-infinite slab. The results obtained 
with this assumption are a good approximation to the solution for the actual body geometry 
when the depth of heat penetration is small compared with pertinent model dimensions, as 
is the case for the present model. A detailed discussion of the accuracy limits due to the 
semi-infinite-slab assumption is presented in reference 27. 
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Test procedure for heat transfer with blowing .- The porous -leading-edge model 
described previously and shown in figure 4 was instrumented to provide heat -transfer 
data on both the porous leading edge and on the thin-skin slab. The locations of the 
thermocouples used in the present investigation are given in table n. 

For the wind-tunnel tests, the desired mass flow rate was established first with the 
model in the retracted position inside the sealed chamber below the tunnel test section. 
(See ref. 20.) During the tunnel startup process this chamber was evacuated to tunnel 
static pressure which was about 0.05 times the pressure on the stagnation line of the 
swept leading edge after the model was exposed to the flow. This sudden change in 
external pressure on the porous leading edge after model injection did not affect the pre- 
set mass flow because of the large pressures (150 to 400 psig) inside the model. 

Heat -transfer data were used during the period from 0.2 to 1.5 seconds after the 
model first entered the air stream. Because of the short testing time, the model was 
practically isothermal; thus heat-conduction effects were minimized and any change in 
the injected mass flow rate which would occur with an increase in model temperature 
(eq. (36)) was minimized. The thermocouple data were recorded at a rate of 20 points 
per second on the digital readout system. 

Data reduction procedures for tests with blowing .- A numerical procedure was 
developed to reduce the temperature -time histories of each thermocouple within the 
porous wall to heating rates. The method provides an approximate solution of the two- 
dimensional heat -conduction problem in a porous matrix with one -dimensional (or radial) 
fluid flow and with temperature -dependent thermal properties. The temperature of the 
fluid is assumed to be the same as that of the porous matrix. Details of this pro- 
cedure and the manner in which it was applied to the porous-leading -edge model are 
given in appendix A. 

Downstream of the interface of the porous leading edge and the slab, the thermo- 
couple wires were spotwelded to the inside surface of the thin stainless -steel skin which 
was 0.030 inch thick. Here, the heating rates were obtained by fitting a second-degree 
curve to the temperature -time data by the method of least squares, over a time interval 
of 1.0 second. The time derivatives of temperature were then computed from the slopes 
of these second-degree curves at the midpoint of the 1.0-second time interval. The 
heating rates, neglecting conduction, were then obtained from the equation 

«w - <p c >m < 27 > 
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RESULTS AND DISCUSSION 


Static Pressure Distribution 

The experimental pressure distributions obtained with the model of figure 2 are 
shown in figure 6. The line faired through the data was used as the distribution of p e 
from which u e and p e were computed for the input to all theoretical solutions 
included herein. An isentropic expansion from p s to p e was assumed for these com- 
putations of u e and p & . 



Figure 6.- Chord-wise distribution of static pressures. 

Zero Blowing 

Heat transfer .- The chordwise distribution of heat -transfer coefficient at the span- 
wise station z/d = 10.6 for Nr 6jCO * 0.92 X 10® is presented in figure 7. The plastic 
models (fig. 3) were used to obtain these data with the phase-change technique. The 
width of the hatched band represents the scatter in 56 data points obtained with eight dif- 
ferent phase-change coatings (0.45 ^ T w /T^ S 0.55) on the three models with different 
thermophysical properties. The theoretical laminar distribution (solid line) was obtained 
with the finite -difference procedure for T w /Tt = 0.4 and e = 0. The good agreement 
between experimental data and the theoretical predictions indicates that the numerical 
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Figure 7-- Comparison of theory and data for laminar flow. 

N Re,oo = 0*92 X 10^; m = 0; z/d = 10.6. 

procedure satisfactorily predicts the experimental heat transfer and that the leading- 
edge boundary layer is laminar at this Reynolds number in agreement with the criteria 
of reference 28. 

Heat -transfer -coefficient contours obtained for typical runs on the plastic models 
with the phase-change coating technique are shown in figures 8(a) and 8(b) for 
Nr 6j oo = 2.6 x 10 5 and 9.3 X 10®, respectively. In figure 8(a) the values of h* decrease 
continuously along the chord at z/d = 10.6. The fact that the contours in this region are 
essentially parallel to the leading edge indicates that the interference disturbances 
between the mounting plate and the model are small; thus, the infinite-cylinder approxi- 
mation is acceptable. 

In contrast to the contours in figure 8(a), the contours at the larger Reynolds num- 
ber (fig. 8(b)) are not parallel to the leading edge at any spanwise station for x/d 1.2. 
On the leading edge for x/d £1.0 and along the chord z/d = 10.6, the contour lines are 
parallel to the stagnation line and the values of h* decrease with increasing x/d. How- 
ever, for x/d 1.3 (at z/d = 10.6) the values of h* increase from a minimum of 
2.8 X 10"® to a maximum of 4.3 X 10"^. At x/d ~ 2.2, the values of h* again begin 
to decrease. The locations of these minimum and maximum contours in h* are too far 
from the mounting plate to be caused by interference effects, as indicated by the results 
of figure 8(a). Oil-flow patterns, to be presented subsequently, substantiate the lack of 
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plate interference effects in this region. Figure 8(a) also shows that the effects of the 
truncated end of the slab are not responsible for the contour behavior noted in figure 8(b) . 

Comparisons and discussion of these different results are best made by means of 
heat -transfer-coefficient distributions. The chordwise distributions of h* along the 
chord z/d = 10.6 for the same test conditions of figure 8 are given in figure 9. Theo- 
retical predictions made with the two eddy- viscosity models described previously (one 
corresponding to the assumptions of invariant turbulence and the other to Lettau's 
(ref. 10) vorticity-transfer hypothesis) for T w /Tt = 0.4 and (Z/6) max and 0.07 

are shown as the dashed lines in figure 9. The solid lines of figure 9 represent predic- 
tions of the laminar theory obtained with e = 0. 

The data and theoretical results in figure 9(a) are for Nr 6)00 = 2.6 x 10®. The 
width of the hatched band represents the scatter in 97 data points obtained with six dif- 
ferent phase-change coatings (0.41 i T w /Tt SI 0.58) on the three models with different 
values of pck. Reference 29 presents comparisons of experimental and theoretical 
heat -transfer results which indicate that the leading-edge boundary layer is turbulent at 
this Reynolds number. In figure 9(a) the two theoretical turbulent distributions (different 
values of (£/6) max ) are within the spread of the experimental data in the leading-edge 
region. The original density of data points for the band shown in this figure indicates 
that the agreement of experimental data and predictions from the finite-difference pro- 
cedure is best when the (Z/ 6) max value of 0.05 is used. On the forward portion of the 
slab (0.785 ^ x/d £ 1.5), the experimental values of h* fall below the theoretical 
turbulent distribution and approach the laminar predictions. The experimental values 
of h* remain slightly above the laminar predictions over the remainder of the slab 
region to x/d = 3.5. 

The experimental and theoretical results for Nr 6>co = 9.3 x 10® are shown in 
figure 9(b). On the leading edge the theoretical turbulent prediction with the invariant- 
turbulence assumption and U/5) max = 0.05 is in best agreement with the data. On the 
forward portion of the slab, the data again drop below the turbulent predictions and tend 
to approach the laminar prediction. However, in contrast with the experimental results 
at the lower Reynolds number (fig. 9(a)), the level of the data begins to increase at 
x/d « 1.3 until agreement with turbulent levels is obtained at x/d « 2.0. 

These trends of the experimental data with x/d and Reynolds number (figs. 9(a) 
and 9(b)) as compared with the theoretical predictions suggest the possibility of laminari- 
zation of the initially turbulent boundary layer on the leading edge. The drop in experi- 
mental heat transfer just aft of the rapid expansion on the leading edge until nearly 
laminar values are reached could also be caused partly by the effect of streamline curva- 
ture on the mixing length. (See refs. 14 and 30.) However, if the effects of curvature 
were the main cause of this drop in heating, it might be expected that the experimental 
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M N Re,~ = 2 * 6 x lo5 « 



0>> %e,« = 9-5 X 1 05. 

Figure 9*- Chordwlse distribution of heat-transfer coefficient, m = 0; z/d = 10.6. 
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levels of heating on the slab would rapidly approach turbulent values at both Reynolds 
numbers since the curvature on the slab is zero. The possibility of laminarization will 
be discussed further after additional data are presented that will confirm the experi- 
mental trends noted in figure 9. 

Surface stre amlines.- Oil-flow patterns obtained on the Slab portion of the pressure 
model at Nj^ e oo = 2.6 x 10^ and 9.3 x 10^ are shown in figures 10(a) and 10(b), respec- 
tively. The streamlines at the wall computed from theoretical values of r w x and r w z 
for both laminar and turbulent flow are included in the figure for comparison. These 
results are referred to in the figures as "surface viscid" streamlines. Inviscid stream- 
lines at the edge of the boundary layer as based on the pressure distribution of figure 10 
and the infinite-cylinder assumption that w e is constant are also shown. The large 
deviations between the viscid and inviscid streamlines indicate that large cross flows are 
present. 

For the lower Reynolds number (fig. 10(a)), the oil-flow patterns are essentially 
parallel to the laminar viscid streamline. This result is in agreement with the heat- 
transfer trends noted in figure 9(a). At Nr 6>00 = 9.3 x 10^ (fig. 10(b)), the oil-flow 
patterns are nearly parallel (except just aft of the leading edge) to the theoretical viscid 
streamline for turbulent flow with the invariant -turbulence eddy-viscosity model. Again 
these results are in agreement with the corresponding trends of figure 9(b). Surface 
oil-flow patterns caused by the boundary-layer pitot -tube rake are evident in figure 10(b). 
Of course, these interference patterns and those near the mounting plate should be dis- 
regarded in the comparisons with theoretical predictions. 

Velocity and Ma ch number profiles. - In order to illustrate the behavior of the veloc- 
ity and stagnation-temperature profiles obtained from the finite-difference solutions, the 

Tf - T 

profile parameters u/u e , w/w e , and — are plotted against y/6 in figure 11. 

' T t,e - T w 

These profiles are in the stagnation region (x/d ~ 0.1) for the two Reynolds numbers of 
2.6 x 10 5 and 9.3 x 10 5 . 

The chordwise velocities "overshoot” the values of u e appreciably for y/6 ^ 0.3. 
These results are qualitatively similar to those of reference 11 where exact solutions for 
the laminar compressible boundary layer on swept cylinders also resulted in overshoot 
of the chordwise velocity profiles. In reference 11 the overshoot phenomenon was 
attributed to the effects of compressibility. Comparison of the computed profiles with 
a power-law profile 
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INVISCID 


■M SURFACE f INVARIANT TURBULENCE ; (Z/S) max =.05 

VISCID V LAMINAR (e = 0 ) 


(a) N Re ,co = 2.6 X 105. 


L- 71-51L 
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^Re,co = 10® 


N Reja) = 2.6X10® 



PROFILE PARAMETER 

Figure 11.- Computed profiles in stagnation region x/d « 0.1. m = 0. 

for n = 1/4 also shown in figure 11 indicates that the present results differ greatly from 
simple power-law profiles, as used for example in reference 3. 

Experimental and theoretical Mach number profiles are shown in figures 12(a) 
and 12(b) for the test Reynolds numbers of 2.6 x 10® and 9.3 x 10®. The experimental 
Mach number profiles were obtained from the ratios of measured pitot pressures and the 
local static pressures given by the faired line in figure 6. The theoretical profiles were 
computed for laminar (e = 0) and turbulent flow. The invariant -turbulence model with 
(Z/5)max = 0-0® was use ^ in the turbulent calculations. 

Although there is considerable scatter in the data, comparisons at the lower 
Reynolds number (fig. 12(a)) indicate good agreement between the experimental results 
and predicted profile shapes and boundary-layer thicknesses for laminar flow. At the 
higher Reynolds number (fig. 12(b)) the experimental results are in excellent agreement 
with the predictions for turbulent flow. The pitot-pressure data are therefore consistent 
with the results of figures 9 and 10 which showed that for x/d > 2.0, the boundary layer 
was nearly laminar at the lower Reynolds number and turbulent at the higher Reynolds 
number. 

Boundary-layer thickness .- The boundary-layer -thickness distributions with x/d 
for Nr 6j0O = 2.6 X 10® and 9.3 x 10® are shown in figures 13(a) and 13(b). The bars 
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(a) \e,«> = 2 ' 6 x 105 * 

Figure 12.- Mach number profiles. All computed profiles for turbulent flow were obtained with the 
invariant-turbulence model and (z/5) max = m = 0. 
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(*) %e,oc = 9.5 X 105. 

Figure 13- - Chordwise distribution of boundary- layer thickness for m = 0 
z/d = 10,6. Invariant -turbulence model used for the two upper lines, 
indicate experimental data spread. 
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indicate the spread in the experimental data. The values of boundary-layer thickness 
were obtained from the Mach number profiles as indicated for typical data in figure 12. 

In figure 13(a) the first data point at x/d = 1.41 is between the lower turbulent predic- 
tion and the laminar prediction. As x/d increases, the experimental values of 5 
approach the laminar predictions up to x/d ~ 3.2, in agreement with the results of the 
heat -transfer and oil-flow data (figs. 9(a) and 10(a)). The last two data points at 
x/d = 4.16 and 5.16 show a slight increase from the theoretical laminar values. Since 
the heat -transfer data of figure 9(a) only extended to x/d = 3.5, any trends for x/d > 3.5 
in the heating data are not available. Close examination of the oil -flow data (fig. 10(a)) 
indicates a slight upward trend near the end of the slab (x/d « 5) toward the turbulent pre- 
diction. However, this trend may be caused partly by interference from the mounting 
plate. 

At the larger value of Nr 6j00 (fig. 13(b)), the first data point for boundary -layer 
thickness (x/d = 1.41) is between the laminar and turbulent distributions. This result is 
again consistent with the corresponding heat-transfer and oil-flow results of figures 9(b) 
and 10(b). The remaining data points are on or slightly above the turbulent result with 
(Z/5) ma x = 0*05 anc * are a ^ so consistent with the heat -transfer and oil-flow data. 

Laminarization .- Comparisons of experimental and theoretical heat -transfer 
results presented herein indicate that the boundary layer on the leading edge of the 
60° swept slab is turbulent for both Nr 6)00 = 2.6 x 10^ and 9.3 x 10*\ These indications 
are consistent with those presented in reference 29. 

The preceding comparisons of data and theory on the slab (x/d > 0.78) for heat 
transfer, oil flow, Mach number profiles, and boundary-layer thickness have all given 
consistent results regarding the apparent trends toward laminar or turbulent flow. To 
provide some further indication of whether these trends were caused by laminarization 
of the initially turbulent boundary layer as the flow expands around the leading edge, 
values of a laminarization parameter have been computed. This parameter is based on 
the form proposed in references 31 and 32 but modified herein for application to three- 
dimensional flow. 


The investigations in references 31 and 32 for conical nozzle flow showed that 
laminarization effects become significant when 


M 


du 


e _e > 2 x io" 6 


Pe u e 


2 dx 


(29) 


This parameter can presumably be extended to three-dimensional flow by using stream- 
line coordinates. That is, by replacing u e with q and dx with ds, the value for 
the parameter of 2 x 10" 6 would be applicable to three-dimensional flow. Thus, when 
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e n e 


(30) 


laminar ization in three-dimensional flows might be expected if the same mechanisms are 
responsible for the phenomenon as in two-dimensional flow. However, calculated values 
of the parameter in relation (30) for the present configuration were at most nearly 
1 order of magnitude smaller than 2 x 10 - ®. Apparently, if the tendency noted in the 
present results for the heat transfer to approach laminar values aft of the leading edge 
were caused by laminarization effects, the criterion of relation (30) would be much 
smaller. A possible explanation of a smaller criterion in three-dimensional flow is as 
follows: It has been established that three-dimensional laminar boundary layers with 
large cross flows are inherently unstable and therefore premature transition due to small 
disturbances is likely. (See ref. 33.) For the present configuration, disturbances gener- 
ated by the leading-edge portion of the mounting plate upstream of the model apex would 
be sufficient to cause transition to turbulent flow at small Reynolds numbers. (See 
refs. 28 and 34.) Consequently, the turbulent boundary layers indicated on the leading 
edge in the present tests may have low levels of turbulence intensity and would therefore 
be highly susceptible to laminarization. The small values of (Z/6) max used in the theory 
(0.05 compared with usual values of approximately 0.09) tend to support the possibility of 
low turbulence levels, since small values of the mixing length cause small values of 
eddy viscosity which imply low levels of turbulence. 


Heat Transfer With Blowing 

Blowing-rate distribu tion on porous leading edge .- The mass-injection-rate dis- 
tribution on the surface of the porous leading edge is given in figure 14 in the form of the 
permeability parameter a. This parameter would be constant if the porous material 
were perfectly uniform in porosity and thickness. The chordwise variation in permea- 
bility is shown in figure 14(a) for the spanwise location z/d = 8.85 and for the range of 
internal pressures shown. The spread of the data at a given x/d location in this fig- 
ure gives an indication of the repeatability of the calibration technique, since the only 
variable was the internal pressure. The chordwise variation of the permeability for a 
range in span of 4 inches (7.85 ^ z/d ^ 11.85) and pressures of 250 and 300 psig is pre- 
sented in figure 14(b). The data spread here is due primarily to the spanwise variation 
in the permeability. The technique used to measure the mass injection rates is described 
in appendix B. 

Application of theory for non similar blowing rates .- The present finite-difference 
procedure can be applied directly to calculation of boundary layers with surface mass 
transfer if the chordwise distribution of the normal velocity v w at the surface is known. 
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(*b) 7.85 ^ z/cL < H. 85 ; p c = 250 and 300 psig. 

Figure 1^-. - Chordwise distributions of permeability on porous leading edge. (Flagged 
symbols indicate data "folded" from negative x/d locations. ) 



The infinite-cylinder condition that spanwise derivatives of all mean-flow quantities are 
8v w 

zero requires that - = 0. This requirement is not met in the present tests because 

of the large spanwise variations in the permeability a as shown by figure 14(b) . Equa- 
tion (36) in appendix B shows that the spanwise variations in v w are directly propor- 

9P W 

tional to the a variations, since — — * 0. 

’ dZ 

Consequently, a direct comparison of theoretical predictions with data is not 
possible. Instead, calculations pertaining to a given value of p c for each of the two 
values of Nj^ e ^ have been carried out for two different blowing-rate distributions 
that represent the maximum and minimum range of a values over the span interval 

7.85 s z/d S 12.85. These results will then be compared with the experimental data over 
the same span interval. Thus the general validity of the theoretical predictions can be 
assessed, provided that the experimental values of heat transfer on both the porous 
leading edge and the downstream slab are affected mainly by leading- edge blowing rates 
within the same span interval. 

The chordwise distributions of the maximum and minimum mass flow measured 
over the span interval 7.85 S z/d ^ 12.85 are shown in figure 15 as the hatched bands 
for values of the internal chamber pressure of 212, 300, and 405 psia. The chordwise 
distributions of v w , corresponding to these maximum and minimum mass flow rates, 
were used as inputs for the computer program. Solutions for these input flow rates were 
then obtained for laminar flow (e = 0) and turbulent flow where the invariant -turbulence 
or the Lettau eddy-viscosity model with (1/6) max = 0.05 was used. Comparisons with 
experimental heat -transfer data for the two test Reynolds numbers of 2.6 x 10® and 
9.3 x 10® for T w /Tt = 0.4 will be presented. The model chamber pressures used for 
these tests were approximately 210 and 300 psia for Nr 6 m = 2,6x 10® and 400 psia 
for N Re,co = 9 - 3x 10°; hence the injection flow rates correspond to those shown in 
figure 15. 

Results for Nj^ e M = 2,6 x 10®.- The predictions of heating rates and the experi- 
mental data for Nj^ e ^ = 2.6 x 10 3 are shown in figures 16(a) and 16(b) for p c = 210 
and 300 psia, respectively. The hatched bands in these figures are bounded by the 
upper and lower limits of the predicted heat-transfer rates for the minimum and max- 
imum limits in m obtained from the boundaries of the shaded bands of figures 15(a) 
and 15(b). Note that the lower limit of m gives the upper limit of heating rate and 
vice versa. The experimental data are shown as the vertical bars. The length of these 
bars represents the variation of heating rates measured within the same span interval 

8.85 S z/d ^ 12.85. Since measurements of recovery temperature with blowing were not 
available, all experimental data and theoretical results with blowing have been normal- 
ized with the local theoretical heating rates for no blowing. These reference values of 
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x/d 

(c) P c = ^05 psia; N Re ^ = 9.3 x 10 5 . 
Figure 15 ♦- Concluded. 
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heating rates for no blowing were computed for the same local values of T w that 
occurred in the test data or were used in the theory with blowing. All reference values 
of heating rate in this and subsequent figures were computed with the invariant - 
turbulence eddy-viscosity model and (Z/5) max = 0.05. 

Comparison of the experimental and theoretical heat-transfer results in fig- 
ure 16(a) indicates that the boundary layer on the leading edge is turbulent and that the 
data agree with predictions for minimum injection rates rather than those for maximum 
injection rates. Just upstream of the slab interface, the experimental heating rates 
decrease and are within the predicted laminar limits. The data on the slab agree with 
laminar predictions for the minimum injection rates. These trends suggest the possi- 
bility of laminarization of the initially turbulent boundary layer on the leading edge, as 
was indicated by the results for no blowing (figs. 9(a), 10(a), and 13(a)). The negative 
heating rates predicted by the theory for laminar flow with maximum blowing rates in the 
vicinity of x/d ~ 0.8 are caused by the large expansion of the flow around the leading 
edge. This large expansion reduces the static temperature of the injected film of air 
near the surface below its temperature farther upstream (which would be almost the wall 
temperature). Hence, when this surface film of air is convected downstream, its tem- 
perature is below the local surface temperature, and heat is then conducted from the 
model to the flow. 

Comparison in figure 16(b) of the experimental and theoretical results for larger 
blowing rates (corresponding to p c ~ 300 psia shown in fig. 15(b)) again indicates that 
the boundary layer on the leading edge was turbulent and the experimental data tend to 
agree with predictions for minimum injection rates. No solution could be obtained for 
the maximum injection rate with laminar flow apparently because "blowoff" or boundary- 
layer separation would have occurred. In contrast with the data of figure 16(a), however, 
the data of figure 16(b) show that turbulent flow persists farther downstream on the slab 
to x/d ~ 2.5. This downstream movement of turbulent flow may be caused by the 
destabilizing influence of the increased blowing rates (corresponding to the increase in 
p c from 210 to 300 psia). 

The rather abrupt drop in heating from x/d = 2.65 to x/d = 2.90 may be caused 
by the smaller upstream blowing rates imposed on the flow over the thermocouples 
located in the range of x/d S 2.90. To clarify this statement, the thermocouple loca- 
tions on the slab are identified in figure 17 according to whether x/d £ 2.90 or 
x/d ^ 2.65. The inviscid and surface viscid (both laminar and turbulent) streamlines 
which leave the porous leading edge at z/d = 7.85 are also shown in the figure. Com- 
parison of these streamline patterns with the two groups of thermocouple locations shows 
that the flow over most of the thermocouple locations for x/d ^ 2.90 originates farther 
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Figure 17.- Thermocouple positions on slab of porous-leading-edge model grouped according 
to whether x/d ^ 2.90 or x/d ^ 2.65 with inviscid and viscid streamlines exiting 
porous leading edge at z/d = 7.85. 

inboard on the leading edge than the flow over the thermocouples for x/d 2 2.65. The 
inboard flow is subjected to smaller blowing rates than the outboard flow, as indicated by 
the data in figure 15(b) for z/d < 7.85. The smaller blowing rates for the inboard flow 
may then account for the trend of reduced heating for x/d £ 2.90. 

The agreement of data on the slab with predicted heat transfer for minimum leading- 
edge blowing is attributed to the reduced levels of injection rates as z/d is reduced. 
(Some of these smaller injection rates for z/d < 7.85 are shown in figs. 15(a) and 15(b).) 
The origin of a given streamline on the porous leading edge is far inboard of the z/d 
position at which it passes onto the slab. (See fig. 17.) The average blowing rate over 
the path of a streamline on the porous leading edge is generally smaller than that mea- 
sured at the outboard z/d stations because of the spanwise variation in permeability. 

The reduction in blowing rates for z/d ^ 7.85 may also affect measured heating 
rates on the leading edge for z/d > 7.85. That is, the local heating rates would tend to 
be larger because of the upstream history of small blowing rates along streamlines. 
Another factor which must be considered is the possible effect of some local blockage of 
the coolant flow in the vicinity of the thermocouple junction. That is, the measured mass 
flow rates at these thermocouple locations tend to be somewhat smaller than most of the 
other measured flow rates within the span interval 7.85 ^ x/d ^ 12.85 as indicated in 
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figure 15(a), where the data points obtained at the junction locations are shown as x 
symbols. Data for the larger blowing rates, to be considered subsequently, indicate that 
local blockage of the coolant flow is not a major factor in the measured heating rates. 

As mentioned previously, the heat -transfer results shown in figure 16(a) for the 
smallest blowing-rate distribution (fig. 15(a)) indicate the same apparent laminarization 
trends noted for the same value of N^ e in the data for no blowing (figs. 9(a), 10(a), 
and 13(a)). Keeping the Reynolds number constant and increasing the blowing rate to that 
shown in figure 15(b) apparently caused turbulent flow to persist farther downstream 
(fig. 16(b)). The reduced heating for x/d > 2.90 for this case may be attributed to the 
history of reduced blowing rates for these streamlines. These results indicate that a 
blowing-rate "threshold" exists. Below this threshold, leading-edge blowing has no sig- 
nificant effect on laminarization. Blowing rates above this threshold may affect the 
response of a turbulent boundary layer to stabilizing influences and delay the develop- 
ment of laminarization trends. 

Results for NR e>00 = 9.3 x 10^.- The chamber pressure for the tests at 

N Re oo = 9.3 x 105 was held at p c ~ 400 psia. The predictions and experimental data 
are shown in figure 18, where the same format for data and theory as that of figure 16 is 



Figure 18.- Effect of blowing on chordwise heat-transfer distribution. Bars 
indicate experimental data spread. = 9*5 X 1C>5; p « 400 psia; 

8.85 ^ z/d ^ 12.85. ’ 


43 


I 


used. Again, no solution could be obtained for the maximum injection rate with laminar 
flow apparently because blowoff or boundary-layer separation would have occurred. The 
predicted ratios of heating rates for the minimum injection flow with the Lettau eddy- 
viscosity model exceed unity over most of the slab because the reference q m _o is f° r 
the invariant -turbulence model. 

Comparison of data and predictions shown in figure 18 again indicates that the flow 
was turbulent on the leading edge with measured heating levels near the predictions for 
minimum blowing. Here the theoretical turbulent prediction with the invariant -turbulence 
assumption and (Z/6) max = 0.05 is in best agreement with the data. The large decrease 
in the measured heating just upstream of the slab is probably again caused by laminariza- 
tion. Following this minimum, the heating on the slab increases rapidly toward the pre- 
dicted turbulent levels for minimum blowing. Both eddy-viscosity models give heating 
predictions within the spread of the data for x/d 2.50. The tendency of the data on the 
slab to agree with predictions for minimum blowing rates at the leading edge is again 
attributed to the reduction in blowing with decreasing z/d. (See fig. 15(c).) These 
results are therefore consistent in every respect with the no-blowing results at the same 
Reynolds number. (Compare figs. 9(b), 10(b), 12(b), and 13(b).) Apparently, this blowing 
level has not significantly affected the laminar ization effects or transition location. The 
tendency of the data on the leading edge to agree with predictions for minimum blowing 
rates is also consistent with the blowing data at the lower Reynolds number (fig. 16) but 
is apparently caused mainly by the spanwise variations in blowing (note data beneath 
shaded bands in fig. 15) rather than by the partial blockage of the injected air at the ther- 
mocouples as suggested by the lower blowing rate distribution (x-symbols in fig. 15(a)), 
since the flow rates at locations of thermocouple junctions (x-symbols in figs. 15(b) 
and 15(c)) indicate practically no partial blockage for these higher blowing distributions. 

Effects of blowing in the stagnation region for turbulent flow.- Solutions with the 
invariant -turbulence eddy-viscosity model with (£/6) max = 0-05 were obtained by the 
finite-difference procedure for several blowing rates in addition to those shown in fig- 
ure 15. The heat -transfer predictions in the vicinity of the stagnation line from these 
solutions are shown in figure 19. The ratios of Ngj-^Ngj- are plotted against the 

blowing parameter B/N st , for which, as in reference 35, the recovery factor with blowing 
for both data and theory was taken as 0.9 of the recovery factor for no blowing. 

The dashed line is a correlation of the computed values to within 1 percent for both 
test Reynolds numbers used in the computations. The experimental data are in reason- 
able agreement with the theoretical prediction over the range of blowing rates available. 
The blowing rates used for the data points were taken as the measured values at the 
corresponding thermocouple locations. 
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Figure 19- - Effect of blowing in the stagnation region 
for turbulent flow, x/d ~ 0.1. 

Also shown for comparison in figure 19 is a correlation in terms of the same 
parameters for flat-plate and cone data. These data are for turbulent boundary layers 
and the correlation is obtained from reference 35. The blowing rates for the present 
data are too small to provide a reliable assessment of the prediction. However, com- 
parison of the correlation lines indicates that blowing near the stagnation line of swept 
cylinders reduces the heating more than a comparable blowing rate on flat plates. 


Cross- Flow Predictions 

As mentioned in the Introduction, the use of an inviscid-streamline coordinate sys- 
tem and the assumption of small cross flow simplify the boundary-layer equations. That 
is, the momentum equation in the streamwise direction and the energy equation become 
independent of the cross flow. If the inviscid streamlines are known, the problem is 
therefore reduced to the solution of an equivalent two-dimensional flow. 

In order for the assumption of small cross flow to be generally valid, the square of 
the ratio of the cross-flow velocity to the streamwise velocity must be small. The maxi- 
mum value of this ratio occurs at the limit as y — 0 (see refs. 5 and 6, for example). 
The small-cross -flow criterion is then 

2 

« 1.0 (31) 

w 

In the present notation this shear-ratio parameter is 

— = tan (0 W - 0 e ) (32) 

P 
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where 8 e is the angle between the spanwise coordinate direction and the local inviscid 

streamline direction, and 0 W = arc tan [-£) . 

\ T z/w 

The chordwise distributions of the small-cross-flow parameter ( T n/ T p^ 2 from 
several of the same solutions used in previous figures are shown in figure 20. The maxi- 
mum values in the curves occur near the interface of the leading edge and slab and range 
from 0.17 to 0.28 for turbulent flow with no blowing. These results are consistent with 
the location and magnitude of maximum values of the cross-flow parameter obtained by 
Bradley (ref. 5) for turbulent flow on a circular cylinder at higher Reynolds numbers. 
These peak values from reference 5 were approximately 0.12. 

The distributions from the solutions for laminar flow with no blowing are essentially 
independent of Reynolds number and the maximum value is 0.44, which is in agreement 
with corresponding values obtained in reference 6. The results given in reference 6 are 
for laminar boundary layers on circular cylinders with the assumption of local similarity. 
Interpolation between the results of reference 6 gives a maximum value for the cross- 
flow parameter of about 0.4 for the present conditions. 

The effect of blowing increases the values of the cross-flow parameter by as much 
as an order of magnitude when the maximum blowing rates are used. For these condi- 
tions, the assumption of small cross flow would not be applicable to the momentum 



Figure 20.- Chordwise distribution of small- cross- flow parameter from theory. 
Invariant- turbulence model with ( l/^)max =0.05 used for turbulent 
computations. 
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characteristics of the boundary layer such as velocity profiles, skin friction, and surface 
streamlines. However, the results of reference 6 indicated that heat-transfer predictions 
from the small cross-flow theory for laminar boundary layers were within 15 percent of 
the correct values when the cross-flow parameter was as large as 3.0. 

CONCLUSIONS 

A finite -difference method has been developed to solve the equations for compress- 
ible turbulent boundary layers on swept infinite cylinders. Predictions of surface heat 
transfer by this method were compared with experimental data on a blunt slab configura- 
tion with and without leading-edge blowing. The leading edge was a hemicylinder which 
was swept 60° with respect to the free-stream flow direction. 

Predictions from the finite-difference method for solving the boundary-layer equa- 
tions were also compared with experimental values of surface streamlines, Mach number 
profiles, boundary-layer thickness, and heat transfer on a similar configuration but with- 
out blowing. The tests were conducted at a free-stream Mach number of 8 and over a 
range of free-stream Reynolds numbers based on leading-edge diameter Nr 6 of 
0.92 X 10 5 to 9.3 X 10 5 . 

Two eddy-viscosity models which utilized general mixing-length functions scaled 
to the boundary-layer thickness were used in the calculations. One was termed the 
"invariant turbulence" model and was a scalar function based on the assumption that eddy 
viscosity is independent of direction. The other model was based on Lettau's vorticity- 
transfer hypothesis (Journal of the Atmospheric Sciences, vol. 21, no. 4, July 1964) and 
was also independent of direction but contained the scalar sum of two velocity-gradient 
terms. 

A numerical procedure for solving the two-dimensional heat conduction problem 
in a porous matrix with one -dimensional fluid flow was also developed. This procedure 
was used to reduce the temperature-time history of the thermocouples near the surface 
of the porous leading edge to surface heat -transfer rates. 

The main conclusions from the investigation are as follows: 

1. The finite-difference method gave reasonable predictions for heat transfer, sur- 
face streamlines, Mach number profiles, and boundary-layer thickness on the test con- 
figuration without blowing for both laminar and turbulent flow. Reasonable predictions 
for heat transfer on both the leading edge and slab were also obtained with nonuniform, 
discontinuous blowing at the leading edge. 
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2. Comparisons of experimental and theoretical heat -transfer results indicate that 
the boundary layer on the leading edge of the 60° swept slab is turbulent for both 
N Re,co = 2.6 x 1q5 and 9 - 3 x 1q5 - 

3. For no blowing, the turbulent boundary layer was apparently laminarized just 
aft of the leading edge. At N^ e ,00 = 2-6 X 10®, the boundary layer remained essentially 
laminar to the end of the slab (5 nose diameters along the chord from the leading edge). 

At Nj^ e rjo = 9.3 x 10® apparent transition back to turbulent flow occurred about 

2 diameters from the leading edge. Predictions from both eddy-viscosity models for 
maximum values of the ratio of the mixing length to the boundary-layer thickness 
(Z/5 )„ qy of 0.05 were then in reasonable agreement with the data. These results were 
based on comparisons of theoretical predictions with all measured quantities mentioned 
above. 

4. For the smallest blowing rates used in this investigation, leading-edge blowing 
had no appreciable effect on either the apparent laminarization phenomena or transition. 
However, for N Re .oo = 2 * 6 x 10®, when the blowing rates were increased by about 40 per- 
cent, turbulent flow apparently persisted to about 2 leading-edge diameters farther down- 
stream than for the small blowing rates. These results indicate that blowing had little 
effect on apparent laminarization until some blowing rate threshold was exceeded. 

5. Comparison of values for a laminarization criterion indicated that the turbulent 
boundary layer on the leading edge for these tests was much more susceptible to laminar- 
ization than a comparable conical nozzle flow. This increased susceptibility to laminari- 
zation is tentatively attributed to the low turbulence levels of the boundary layer on the 
leading edge as evidenced by the relatively small values of (Z/6) max required to obtain 
agreement of theory with experiment. 

6. Measured heat -transfer rates on the leading edge over a 4-inch span interval 
agreed with values predicted for the minimum injection rates measured on the leading 
edge over the same 4-inch interval. This agreement with predictions for minimum 
injection was attributed mainly to the smaller injection rates near the apex of the leading 
edge. The influence of the partial blockage of injection at the thermocouple locations on 
the leading edge was considered minor. 

7. Downstream of the leading edge on the thin-skinned slab, the measured heating 
rates were also generally in agreement with predictions for the minimum injection rates 
on the leading edge. These results are attributed to the smaller injection rates near the 
apex of the leading edge due to the spanwise distribution of permeability. 
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8. Blowing on the leading edge increased the magnitude of the cross-flow velocities 
by as much as an order of magnitude. Therefore, the assumption of small cross flow, 
which was not used in the present method, would not be expected to provide reliable pre- 
dictions for the momentum characteristics of boundary layers with large leading-edge 
blowing. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 28, 1971. 
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APPENDIX A 


A HEAT-BALANCE METHOD FOR THE INVERSE SOLUTION OF 
TWO-DIMENSIONAL HEAT CONDUCTION IN A POROUS 
MATRIX WITH ONE -DIMENSIONAL FLUID FLOW 

A numerical procedure was developed to reduce the temperature-time histories of 
each thermocouple within the porous wall to heating rates. The method provides an 
approximate solution of the two-dimensional heat -conduction problem in a porous matrix 
with one -dimensional (or radial) fluid flow and with temperature-dependent thermal 
properties. The temperature of the fluid is assumed to be the same as that of the porous 
matrix. 

To illustrate the procedure, consider the finite slab shown divided into six blocks in 
figure 21. The upper and lower surfaces of the slab are insulated. The initial tempera- 
ture and time -dependent one -dimensional blowing rates for each block are known. At time 
t a t 0 , the front face of the slab (the left side in the figure) is exposed to an unknown heating 
rate which may vary with time and block location. The temperature change with time at 



BLOCK CONFIGURATION WITH HEAT AND MASS FLUX (^=0) 



TEMPERATURE DISTRIBUTION 


Figure 21.- Sketch of six-block heat-conduction model to illustrate method and 
notation. Model is assumed two dimensional; that is, derivatives of all 
quantities in z-direction are zero. 
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APPENDIX A - Continued 


the midpoint of both surface blocks is measured. In order to determine the temperature - 
time profiles of the remaining blocks, a heat balance is written for each. In the heat 
stored term each block is assumed to be at a uniform temperature throughout its volume. 
In the heat conduction and gas absorption terms, a linear temperature variation between 
centers of blocks is assumed. The heat balance for block 1 is as follows: 


Con- 

vection 


*1*4,1 - \ 


Radiation 

W'tv ) 4 

x conducted, gas 
2 (%)l,2 A x, 


. x conducted, 

Stored Gas absorbed material 

(T, - T'A x.(Tj - T„) . (T. - T 2 ) 

Pl c l V l At ' ‘ 2Al 3 A x,3 C p,l (Xj + x 2 ) " 2k l,2 A x,l ( Xl + x 2 ) 

y conducted, material y conducted, gas 
2 (kg)l,4 A y,l \y\ + y A ( " 2 (^s)l, 




(33) 


where the primes denote a temperature evaluated at the previous time step and the sub- 
scripts 1 to 6 refer to block numbers. The subscript q indicates the surface exposed 
to heating rate q. Since the calculation starts at t = t Q , all the primed quantities are 
known. A similar heat balance is written for all blocks and the resulting system of 
linear equations are 


a l T l + b l T 2 


+ d l T 4 


= “ A q,l 4 l + g l 




a 2 T l + b 2 T 2 + c 2 T 3 


b 3 T 2 + c 3 T 3 


a 4 T l 


b 5 T 2 


+ c 6 T 3 


+ e 2 T 5 

— 

g 2 


+ f 3 T 6 = 

g 3 

+ d 5 T 4 + e 4 T 5 

= 

g 4 

+ d 5 T 4 + e 5 T 5 

+ f 5 T 6 * 

g 5 

+ e 6 T 5 

+ f 6 T 6 - 

g 6 


(34) 


J 


The temperatures T-^ and are known; that is, they are obtained from the mea- 

sured temperatures. Also, the temperature coefficients a^, b^, . . f^ and the con- 
stants g^ are calculated from the given inputs. The unknowns are then q^, T 2 , 

Tg, Tg, and Tg which may be obtained from the solution of the following determinant 
system: 
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Unknowns 

^1 t 2 T 3 ^4 T 5 T 6 
Determinant elements 


A q,i b i ° 

0 b 2 c 2 

0 b 3 C 3 

0 0 0 

0 b 5 0 

0 0 c 6 


0 0 0 

0 e 2 o 

0 0 f 3 

A 4,4 e 4 0 

0 e 5 0 

0 e 6 f 6 




Constants for any n time step 
_a l T l " d l T 4 + g l 
~ a 2 T l + g 2 

+ g 3 

" a 4 T l " d 4 T 4 + g 4 
- d 5 T 4 +g 5 


( 35 ) 


This system is solved successively at each n time step where t = t D + nAt. The values 
of the specific heats and thermal conductivities are "updated" at the new temperatures 
obtained after each time step. The coefficients and constants of the determinant can then 
be updated, and the heating rates and temperature distributions obtained for the next time 
step. This procedure is repeated until a specified time interval tf = t D + jAt is reached. 
The surface heating rates for each block are thereby computed over this specified time 
interval with the time -dependent inputs Tj., T 4 , mg, and nig and the temperature- 
dependent material properties. 

Two different block configurations were used for five spanwise stations on the por- 
ous leading edge, as shown in figure 22. These block configurations were set up so that 
experimental thermocouple readings could be used as direct inputs to the computer pro- 
gram without smoothing techniques applied to the data. Hence, the chordwise thermo- 
couple distributions at the spanwise stations as indicated in the figure determined the 
block configurations because the temperature-time history for one block in each column 
of blocks is required as an input to the program. In the 56-block configuration, the sym- 
metry of the thermocouple locations about the x/d = 0 station allowed the use of the 
"surface" thermocouple readings at both the location of the thermocouple and its "mirror" 
image. This increase in the number of columns from four to seven was made in order to 
reduce the chordwise length of the blocks and still provide thermocouple data in the chord- 
wise center of each column of blocks. Superposition of the thermocouple locations in this 
manner is strictly valid only if the chordwise distribution of permeability a is perfectly 
symmetrical. Deviations in a from symmetry can be obtained from figure 14 where the 
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o* 



0 ° 



Figure 22.- Block configurations used to reduce data from porous leading 
edge. Thermocouple locations are given in table II and designated as 
x in sketch. 

data for the x locations of 10°, 30°, 40°, and 70° are shown as flagged symbols. These 
deviations are small for most of the data and hence were neglected in the 56-block 
configuration. 

The chordwise conduction errors are small, in spite of the relatively large block 
dimensions in the chordwise plane, because of the nearly isothermal conditions during 
the short test period. The boundary condition of the insulated edge (fig. 22) will cause a 
slight reduction in the heating rates calculated for the two edge blocks compared with that 
which occurred in the tests since no such restriction applied to the experimental model. 

The surface block of each column was made 0.004 inch thick in accordance with the 
approximate depth of the thermocouple junction. The remaining thickness of approxi- 
mately 0.090 inch was divided into seven blocks of equal thickness. The time step used 
herein for all solutions was 0.05 second. Comparisons of results from calculations for 
one-dimensional heat-conduction problems with exact solutions from Carslaw and Jaeger 
(ref. 36) indicated that for the block thicknesses and time steps of the present computa- 
tions, errors would have been negligible had only five equal-thickness blocks been used. 
To assess further the reliability of the heat-balance method, tests were conducted on the 
porous -leading-edge model with no mass injection at the same tunnel conditions for which 
the phase-change heat -transfer technique was used on the plastic models. The validity 
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of the inverse heat-balance method was confirmed by comparisons of the resultant heat- 
transfer distributions. Description and verification of a one -dimensional inverse heat- 
balance solution without fluid flow is given in reference 37. 
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TECHNIQUE FOR MEASURING BLOWING-RATE DISTRIBUTION 
ON POROUS LEADING EDGE 

The local mass injection rates on the porous leading edge were measured at several 
spanwise stations and chordwise distances over a range of internal pressures from 150 
to 400 psig. These measurements were made with a plastic tip probe (0.096 inch inside 
diameter) which was pressed against the porous surface as shown in figure 23. The probe 



Figure 2J.- Plastic tip probe and alinement jig for measuring injection rates. 

was connected with a butyl phthalate manometer and a pressure-sealed container of known 
volume. The total mass flow rate through the porous leading edge was set at the desired 
level by adjusting the internal pressure. This flow rate was measured with a floating 
ball-type flowmeter. The plastic tip probe was held in contact with the model for a mea- 
sured period of time. At the end of this period (§30 sec) the system was sealed just 
downstream of the plastic tip probe. After equilibrium was established, the change in 
pressure and volume (the change in volume due to the change in level of the manometer 
fluid) of the isothermal system was recorded. The local mass flow rate through the con- 
tact area of 0.096-inch diameter could then be calculated. The pressure rise in the 
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sealed volume was limited to less than 0.7 percent of the pressure inside the model. 
Therefore, the mass-flow distribution through the porous material was not altered by the 
measurement technique. The validity of this statement is evident from the relation for 
the mass flow rate through a uniformly porous material. This relation is given in refer- 
ence 38 (eq. (7), p. 78) and, in the notation herein, is 


m = P^ v w = 


a 


w 


w w 2R(pT) w AY 


(36) 


Measurements were made at chordwise intervals of 10° around the leading edge 
and at spanwise intervals of 1 inch or less along the porous leading edge in the region 
7.85 ^ z/d S 11.85. 
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